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Abstract 

While the ecology and evolution of partial migratory systems (defined broadly to 
include skip spawning) have been well studied, we are only beginning to under- 
stand how partial migratory populations are responding to ongoing environmen- 
tal change. Environmental change can lead to differences in the fitness of 
residents and migrants, which could eventually lead to changes in the frequency 
of the strategies in the overall population. Here, we address questions concerning 
the life history of the endangered Gila cypha (humpback chub) in the regulated 
Colorado River and the unregulated tributary and primary spawning area, the 
Little Colorado River. We develop eight multistate models for the population 
based on three movement hypotheses, in which states are defined in terms of fish 
size classes and river locations. We fit these models to mark-recapture data col- 
lected in 2009-2012. We compare survival and growth estimates between the Col- 
orado River and Little Colorado River and calculate abundances for all size 
classes. The best model supports the hypotheses that larger adults spawn more 
frequently than smaller adults, that there are residents in the spawning grounds, 
and that juveniles move out of the Little Colorado River in large numbers during 
the monsoon season (July-September). Monthly survival rates for G. cypha in 
the Colorado River are higher than in the Little Colorado River in all size classes; 
however, growth is slower. While the hypothetical life histories of life-long resi- 
dents in the Little Colorado River and partial migrants spending most of its time 
in the Colorado River are very different, they lead to roughly similar fitness 
expectations when we used expected number of spawns as a proxy. However, 
more research is needed because our study period covers a period of years when 
conditions in the Colorado River for G. cypha are likely to have been better than 
has been typical over the last few decades. 



Introduction 

Partial migration, in which a portion of a population is 
migratory and the rest is sedentary, is a widespread phe- 
nomenon observed in a variety of taxa, including fish 
(Jonsson and Jonsson 1993), birds (Lundberg 1988), 
amphibians (Grayson et al. 2011), and reptiles (Blake 
et al. 2013). Partial migratory systems are complex and 
may include individuals with movement strategies that 
are not easily classified. Nonetheless, three idealized forms 
of partial migration are recognized: (1) shared breeding 



(spawning) areas: residents and migrants occupy the same 
general area during breeding season and migrants move 
away from the breeding grounds during the nonbreeding 
season (e.g.. Lack 1943), (2) shared nonbreeding areas: 
residents and migrants occupy the same general area dur- 
ing the nonbreeding season and migrants move elsewhere 
to breed (e.g., Morrissey et al. 2004), and (3) skip spawn- 
ing: only a fraction of potential breeders leave nonbreed- 
ing areas in a given year to travel to breeding (or 
spawning) grounds (Shaw and Levin 2011). Although par- 
tial migration of all three types is widespread, little is 
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known about how human modification of the environ- 
ment is impacting different portions of partial migratory 
populations. Existing studies suggest that impacts of envi- 
ronmental change on partial migratory systems can be 
nonintuitive (Berthold 2001; Nilsson et al. 2006) and 
involve both positive and negative impacts on the fitness 
of a given strategy (Hebblewhite and Merrill 2007, 2011). 

River systems are among the most modified ecosystems 
worldwide (Poff and Zimmerman 2010), and dams can 
have profound impacts on fully and partially migratory fish 
populations (Pringle 2003). Much research has focused on 
impacts dams have on connectivity in river systems such as 
severing some migration routes and making others more 
arduous (Ward and Standford 1995; Pringle 2003). In 
many cases, dams also change downriver environmental 
conditions (e.g., temperature; Olden and Naiman 2010) 
and are associated with changes in ecological communities, 
including introductions of nonnative fish (Ward and 
Standford 1983; Cross et al. 2013). Changes in environ- 
mental conditions impact vital rates for populations that 
spend some or all of their time in these waters (Poff and 
Hart 2002). For example, Nelson et al. (2002) concluded 
that the decline of the migratory contingent of partial 
migratory populations of Salvelinus confluentus (buU trout) 
found in the headwaters of the Bitterroot River (Montana) 
was not directly related to presence of the dam, but rather 
to changes in factors such as predation or temperature in 
the main river (hereafter mainstem). 

Prior to introduction of nonnative fishes and completion 
of numerous dams, federally listed Gila cypha (humpback 
chub, hereafter chub) was widely distributed in the Colo- 
rado River (USFWS 1994). Chub have since declined 
throughout their range and are restricted to six extant pop- 
ulations. The largest of these populations, and the focus of 
our research, is isolated from upstream populations by Glen 
Canyon Dam. Over 95% of this population resides in the 
lower 14 km of the Little Colorado River (hereafter LCR) 
and in the mainstem Colorado River within a few km of the 
LCR (Valdez and Ryel 1995; Fig. 1). Chub population 
trends are often attributed to conditions in the Colorado 
River portion of this river system, which is regulated and 
highly altered in terms of temperature, turbidity, and flow 
and which contains many introduced species, including sal- 
monids. This population declined during the late 1990s, 
coincident with cooler water temperatures and higher sal- 
monid abundances, and partially recovered over the last 
decade, when water temperatures were warmer and salmo- 
nid abundances were lower (Coggins et al. 2006a,b; Yard 
et al. 2011). 

Chub spawn almost exclusively in the LCR during April 
and May. Juvenile chub that move from the LCR to the 
Colorado River at small sizes in May or June rarely sur- 
vive (Robinson et al. 1998; Robinson and Childs 2001); 
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however, the fate of juveniles leaving the LCR during the 
monsoon season between July and September appears to 
vary with environmental conditions (e.g., water tempera- 
ture and salmonid abundances). Research during the 
1990s suggested that the number of juveniles leaving the 
LCR between July and September can be substantial; 
however, very few of the fish appeared to survive (Valdez 
and Ryel 1995; Clarkson and Childs 2000). Given the 
presence of adult chub in the Colorado River, many sci- 
entists and managers in the system concluded that suc- 
cessful chub likely reared in the LCR before switching to 
a migratory life history at larger sizes. More recent studies 
suggest that juvenile chub can rear in both the LCR and 
Colorado River (Finch 2012), resurrecting basic questions 
about chub life history. The rate of this juvenile outmi- 
gration is unknown, as are the rates of movements 
between the Colorado River and LCR at later life stages. 
Adults that rear in the Colorado River eventually return 
to spawn; however, it is believed that some unknown pro- 
portion of the adult population forgoes spawning every 
year (Coggins et al. 2006b). In addition, Douglas and 
Marsh (1996) showed that adults can be found in the 
LCR in all months of the year and hypothesized that a 
portion of the adult population might reside in the LCR 
year-round. That is, the population may exhibit both a 
shared breeding partial migratory system and a skip- 
spawning system. 

Here, we ask four general questions about chub life his- 
tory. Firstly, what is the rate of skip spawning in chub 
and does this rate change as fish grow? A priori, we 
hypothesized that larger adults should spawn more fre- 
quently, as has been suggested for chub (Gorman and 
Stone 1999) and observed in other species (Secor 2007). 
Secondly, are there chub that reside in the LCR year- 
round? Thirdly, what are the rates of movement between 
the LCR and the Colorado River prior to adulthood? 
Lastly, how do survival and growth of different size clas- 
ses vary between the Colorado River and the LCR and 
what is the fitness of migrants relative to residents? Previ- 
ous work suggests that fish grow more quickly in the 
LCR; however, comparisons of survival rates for different 
portions of the river have not previously been estimable 
because of a lack of systematic sampling in the Colorado 
River. We also estimate population size for different size 
classes of chub, including the total adult abundance and 
link trends in abundance to movement hypotheses. 

Materials and Methods 
study sites 

The chub population that spawns in the Little Colorado 
River is primarily concentrated in the 13.6 km of the 
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LCR that are closest to the Colorado River, and an adja- 
cent ~17 km section of the Colorado River referred to as 
the LCR aggregation (Valdez and Ryel 1995). The Colo- 
rado River is colder (average difference of 6.4°C during 
the study period) and -50 times larger than the LCR (see 
uses gauges 09402500 and 09402300 for detailed dis- 
charge and temperature records). Our analysis is based on 
data collected in the LCR study site, which spans the 
lower 13.6 km of the LCR, and the Colorado River study 
site, which spans ~2.7 km of the Colorado River just 
downriver of the confluence of the LCR and Colorado 
Rivers (Fig. 1 & Table 1). 

Data collection 

Fish in the Colorado River study site were collected using 
two gear types: 47-90 un-baited hoop-nets (50-60 cm 
diameter, 100 cm long, single 10-cm throat, made of 6- 
mm nylon mesh, checked daily for 8-12 days per trip) 
and slow-speed boat electrofishing (pulsed DC current, 
15-20 amps, 200-300 volts, boat speed 7-10 sec per 
meter of shoreline, repeated 24-72 h apart for 3-5 total 
passes per trip). Sampling in the LCR consisted of mark- 
recapture events spanning the whole LCR study site and 
mark-only events over a more limited spatial extent 
(Table 1). Mark-recapture events in the LCR included 



three passes at 180 net locations using un-baited hoop- 
nets (for more details see Van Haverbeke et al. 2013). 
Mark-only events occurred between 0.8 and 1.7 km upri- 
ver of the confluence in 2009-2011 and between 1.8 and 
3 km in 2012. In 2009-2011, fish were caught using 7-10 
hoop-nets over the course of a few days. In 2012, sam- 
pling occurred over 4 days using 18 hoop-nets and was 
augmented by targeted sampling for juveniles using seines 
and dip nets. 

Upon capture, chub were measured, scanned for prior 
marks, and, if appropriate, marked. All chub >100 mm 
total length (TL) that did not already have a 134.2 kHz 
passive integrated transponder (PIT) tag received a PIT 
tag with a unique number allowing future identification 
of individual fish. Juvenile chub (40-100 mm TL) 
received a visible implant elastomer (VIE) mark that 
identified the current sampling trip and spatial location. 
This occurred regardless of whether a fish was previously 
marked, so juveniles accumulate marks if they are caught 
on multiple trips until they grow to greater than 100 mm 
TL. As a consequence, each juvenile carries its capture 
history with it (in the form of VIE marks); however, 
information on individual covariates is lost. Juveniles that 
were captured or recaptured in the LCR study site during 
spring sampling received no mark, and these captures 
were not included in our analysis. 
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Table 1. Approximate timing of important biological events and sampling efforts. 



Month 



J 



Biological events 
Migrating chub move from Colorado River to LCR 
Chub spawn 

Most adult chub leave LCR and return to Colorado River 
Monsoon season 
Sampling efforts 
System-wide mark-recapture in LCR (2009-2012) 
Spatially restricted mark-only in LCR (2009-2011) 
Spatially restricted mark-only in LCR (2012) 
Sampling in Colorado River study site (2009-201 1) 
Sampling in Colorado River study site (2012) 



Grey shading indicates the approximate timing of biological events and X's indicate months during which sampling occurred. 



Model overview 

We use multistate models conditioned on first capture 
(Arnason 1973; Brownie et al. 1993; Schwarz et al. 1993; 
Lebreton et al. 2009) to test movement hypotheses and 
determine the relative fitness of different life-history 
strategies. Multistate models allow ecologists to estimate 
transitions between biological states (e.g., size or age 
classes, spatial location) when capture and/or survival 
probabilities vary between states and are less than one. 
In our models, states are defined in terms of the size of 
chub and their capture location and we estimate move- 
ment and growth rates, in addition to survival and cap- 
ture probabilities (Fig. 2). We define five size classes and 
test hypotheses about the movement of juveniles (size 
class 1: 40-100 mm TL), smaller adults (size class 4: 
200-250 mm), and larger adults (size class 5: >250 mm). 
We consider three locations: (1) the LCR study site, (2) 
the Colorado River study site (Fig. 1), and (3) the Colo- 
rado River outside of the monitoring site. States one 
through five correspond to the five size classes in the 
LCR, states six through ten correspond to the five size 
classes in the Colorado River study site, and states 11 
through 15 correspond to the five size classes elsewhere 
in the Colorado River (Fig. 2). Inclusion of states 11 
through 15 allows us to account for fish in the Colorado 
River outside of the Colorado River study site that tem- 
porally emigrate into the LCR to spawn some years, but 
are otherwise unavailable for capture. Chub vital rates 
(growth, survival, etc.) for a given size class in the Colo- 
rado River study site are assumed to be representative of 
vital rates for all chub in that size class in the Colorado 
River (Fig. 2). 

Our study period is from April 2009 and October 
2012, and we assume a monthly time step for modeling. 
Of a possible 43 months, the Colorado River study site is 
sampled 15 times and the LCR is sampled 23 times 



(Table 1). For computational reasons (our data include 
30,067 individual capture histories) and because we did 
not seek to fit models with between-trip trap responses or 
individual variation, we decomposed capture histories 
into a series of release-recapture events and determined 
the frequency of these events (sensu Pradel 1993). In 
other words, we summarized the information in our indi- 
vidual capture histories in an m|J array (Williams et al. 
2002, p. 456) and analyzed the information contained in 
this array. 

Survival and state dynamics 

In our model survival, (f>\ is defined as the probability 
that a marked fish in state I at time t survives from 
month f to f -L 1, regardless of whether the individual 
transitions into another state during the same interval. 
All the models examined in this paper assume that sur- 
vival is constant through time (i.e., = </>'), but allow 
survival to vary between states (with the exception that 
fish in the same size class in the Colorado River have 
the same survival regardless of whether they are in the 
Colorado River study site or elsewhere). Our code is 
written to allow future examination of variation in sur- 
vival over time (Appendix SI). State dynamics (i.e., 
movement and growth) were characterized generally 
using state transition parameters, (/'[''", reflecting the 
probability that a fish in state I at time t will transition 
to state m at time t + I, conditional on survival until 
t + I. For the purposes of our study, transition probabil- 
ities are expressed in terms of (1) growth conditional on 
survival and (2) movement conditional on survival and 
grovrth. Growth was assumed to be constant over time, 
to vary between those states in which growth was possi- 
ble (i.e., it is not possible to transition out of states 
associated with size class 5), and to be the same for a 
given size class found in either of the two Colorado 
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Colorado River 
study site 



Little Colorado River study site 
(Tributary - Spawning) 



Rest of the Colorado 
River (Unobserved) 



Juvenile 
Size class 1 
(40 - 100 mm) 



Small sub-adult 

Size class 2 
(100 - 150 mm) 



Large sub -adult 

Size class 3 
(150 - 200 mm) 



Small adult 
Size class 4 
(200 - 250 mm) 




Large adult 
Size class 5 
(250 mm+) 

) / — ' — *• 

Summary of model parameters 

(f) - survival O) - movement y - size transition (growth) 

T - proportion of Colorado River fish in Colorado River HBC monitoring site 

Figure 2. Schematic representation of the general form of the multistate model used in this paper. The number on each fish is its state. 



River locations. As such, we considered eight growth 
probabilities, where the first four parameters 

= 1,2,3, or 4) correspond to the probability of transi- 
tioning out of size classes one through four in the LCR, 
and the next four (z = 5,6,7, or 8) are defined as the 
probabilities of growing out of size classes one through 
four in the Colorado River. 

Movement probabilities were constrained based on 
two assumptions. First, we assumed that movement of 
chub between the Colorado River study site and other 
parts of the Colorado River is negligible. We base this 
assumption on studies over both short and long tempo- 
ral scales that suggest that chub exhibit strong patterns 
of spatial fidelity (outside of migration) when they are 
in the Colorado River, restricting their movements to an 
average linear extent of around 500 m (Valdez and Ryel 
1995; Paukert et al. 2006; Gerig 2012). Second, we 



assumed that a constant proportion, t, of chub within 
the Colorado River are found within the Colorado River 
study site, regardless of size class. This assumption 
should be approximately true provided that vital rates in 
the Colorado River study site are similar to elsewhere in 
the Colorado River. Based on these assumptions, we 
decomposed movement in a vector or probabilities, 
representing the probability of a state transition corre- 
sponding to moving into the same size class in the other 
part of the river system (i.e., moving in the Colorado 
River if currently in the LCR or vice-versa) between time 
t and f -h 1. For chub moving from the LCR into the 
Colorado River, coj was multiplied by t to express move- 
ment into the Colorado River study site and multiplied 
by (1 — t) for movement into areas outside of the Col- 
orado River study site. For states associated with size 
classes 2 and 3 (i.e., 1 = 2, 3, 7, 8, 12, or 13), a)\ was 
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assumed to be constant through time; however, var- 
ied by month (and depending on prior capture history) 
for other size classes based on a priori hypotheses (see 
below). As with growth and survival, movement proba- 
bilities out of the Colorado River were assumed to be 
the same regardless of whether a fish was in the Colo- 
rado River study site or elsewhere in the Colorado (e.g., 
0/ = «'^). Lastly, because no instances of juvenile fish 
moving from the Colorado River to the LCR were 
observed in our data set, the associated parameters (co*) 
were fixed at zero, for all f. 

We incorporated the joint effects of survival and transi- 
tion probabilities in an array, (^J'", to express the proba- 
bility that an individual in state I at time t both survives 
and transitions to state m in time f -I- 1. For example, the 
probability of moving from state 1 (juvenile in LCR) to 
state 6 (juvenile in the Colorado River study site) can be 
expressed as follows: 

c;'^ = 0Sa;''' = 0'(i-7')«;t (d 

which states that the probability of surviving and moving 
from state 1 to state 6 is the product of the probability of 
surviving in state 1 (0'), the probability of not growing 
into state 2 (!—"/), the probability of moving from the 
LCR into the Colorado River (coj), and the probability of 
ending in the Colorado River study site (t). 

Capture probability and tag loss 

We assume that when states were observed, they were 
observed perfectly (i.e., fish were assigned to the correct 
size class). We account for the possibility that chub could 
be in a particular state and be unobserved using a vector of 
parameters, p', indicating the probability that a fish in true 
state / during time t is captured. AU of the models we fit 
estimate capture probability separately for each state and 
time period in which there was a possibility of recapture. 
As states 11 through 15 are never directly observed, capture 
probabilities for these states are always set to zero. To 
account for tagging mortality and tag-shedding rates, we 
incorporated additional parameters into our models and 
fixed them at values estimated from laboratory experi- 
ments on sister taxa [Gila elegans and Gila intermedia; 
Appendix S2; Ward et al. 2008). For both PIT and VIE 
marks, the instantaneous rate, d, of tag mortality and tag 
shedding was estimated at 3%. This rate represents mortal- 
ity and tag loss within a few days of tagging and was only 
applied once for newly tagged fish. In addition, as our lab- 
oratory studies (Appendix S2) show that VIE loss rate is 
higher for individuals growing quickly in warmer water, 
we also include a monthly mark loss rate (71) of 0.03 for 
VIE tagged fish in the LCR where juvenile growth rates are 
higher. 



A priori hypotheses 

We evaluate a set of eight models to address three move- 
ment hypotheses (Table 2). Hypothesis 1 (HI): Larger 
adults spawn more frequently than smaller adults, so 
chub in state nine (size class four) will have lower move- 
ment rates (ffl^^ ^ where the subscript refers to the inter- 
val between March and April) from the Colorado River 
to the LCR between March and April than chub in state 
ten (cuj^). H2: Adult chub caught in the LCR during the 
fall may be residents [R) and thus more likely to remain 
in the LCR in the future than other adults (likely 
migrants, M). In other words, movement of potential 
adult residents out of the LCR (oj** and co^^ if rates are 
estimated separately for smaller and larger adults) will be 
less than movements of migrants, especially in the 
2 months after spawning (cJ^ and - where the sub- 
scripts refer to the period between May and June). H3: 
There is a large outmigration of juvenile fish from the 
LCR during the monsoon season in their first year of life, 
so models that estimate a different rate of movement out 
of the LCR for juvenile chub between July and September, 
fflfg, than during the rest of the year, ojj,, will fit the data 
better than a model that assumes a constant movement 
rate, co', throughout the year. Although we considered 
models that allowed for variation between months 
according to our a priori hypotheses, all models assumed 
that there was no interannual variation in rates. 

In addition to addressing the three movement hypothe- 
ses via model comparison, we are also interested in the 
degree to which estimates of chub survival and growth 
for different size classes vary between the Colorado River 
and the LCR. A priori, we hypothesized that growth rates 
would be lower in the Colorado River than in the LCR 
when comparing within a size class and that growth rates 
would decline for larger size classes. We also hypothesized 
that survival would increase for larger size classes, but 
had no a priori hypotheses about how survival would 
compare between the Colorado River and the LCR. Lastly, 
we hypothesized that t would have a value of 0.15-0.25 
based on past studies suggesting that most chub in the 
Colorado River are found near to the LCR confluence. 

Models are fitted by maximum likelihood using R. 
Overdispersion was estimated using the c statistic com- 
puted using a scaled Pearson's chi-squared test statistic 
(Burnham and Anderson 1998, page 68). We compared 
models using quasi-AIC and used c to adjust variance 
estimates for parameters. 

Derived parameters 

To illustrate the consequences of the differences in survival 
and growth rates between the Colorado River and LCR, we 
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Estimated movement parameters that differed between models 
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caught in the 


movement out 






Model 


(in addition to to^, co^, o7, co^, which were 


smaller and 


ICR in fall (i.e., 


of ICR during 






Number 


estimated in all models)' 


larger adults? 


ICR residency)? 


monsoon season? 




AQAIC^ 


1 


co],. col; <r. »r: coISa. 


Y 


Y 


Y 


178 


0 


2 


a^t: o,'"; co»f,: o/,"; c«^7; o.™ 
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N 
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58 


3 




N 
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Y 


173 


91 
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Y 
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5 
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N 
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153 
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Y 
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168 


7 
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CO , CU^ft, Oil, O*.. fOMA. f". . O^MI. 


Y 


N 


N 


175 


186 


8 




N 


N 


N 
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228 



'The subscript for each to refers to the between month intervals for which the parameter was estimated and the superscript refers to the state 
from which movement occurred. Codes in the subscript are as follows: MA - March to April; MJ - May to July; JS - July to September. A w with- 
out a subscript indicates that the parameters was constant across all months, while an * in the subscript indicates that a parameter was estimated 
for all other intervals besides the interval(s) with a special parameter for that state. Superscripts that include two states separated by an "&" indi- 
cate that this parameter is restricted to be the same for both of the states. Superscripts including a letter refer to a separate parameter being esti- 
mated for fish with ("R" - resident) and without ("M" - migrant) a prior capture during a fall LCR sampling event. 

is the total number of estimated parameters in the model. 
^AQAIC is the difference in quasi-likelihood between each model and the best overall model in the model set. 



calculate a number of derived parameters based on the 
maximum-likelihood estimates of survival and size transi- 
tion rates and assuming that fish strictly adhere to either a 
skip-spawning migratory or breeding area resident strategy. 
In presenting these analyses, we acknowledge that some 
fish do not fall neatly into either category, that the period 
of our monitoring almost certainly does not reflect the 
likely variation in vital rates, and that we are ignoring vari- 
ous sources of uncertainty. The specific parameters we 
derive are the expected number of months spent in each 
size class, the probability of surviving from July of a fish's 
first year to the smaller mature size class (size class 4), the 
probability of surviving from July of a fish's first year to 
the larger mature size class (size class 5), and the expected 
life span of size class 5 chub conditional on them having 
survived to that size (and in the case of migratory fish, 
depending on whether they spawn every year, skip spawn, 
or never spawn). As a measure of the relative fitness of 
each strategy, we also derive the expected number of 
spawns as a large adult for a fish that survives to July of its 
first year and then adopts either the migratory or the resi- 
dent strategy by multiplying the probability of surviving to 
size class 5 under each strategy, by the expected number of 
spawns after reach size class 5. We also derived estimates of 
abundance for different states and for the overall adult 
abundance using the Horvitz-Thompson estimator (see 
Appendix S3 for details). 



Results 

The best model according to QAIC supported all three 
movement hypotheses (Table 2) and had a c of 1.7, indi- 
cating slight overdispersion (Appendix S4). The best model 
suggests that some adults skip spawn and larger adults skip 
less fi-equently. In others words, larger adults have a much 
greater probability than smaller adults of moving into the 
LCR from the Colorado River in the month prior to 
spawning (0.61 vs. 0.31; Fig 3A). In addition, adults of 
both size classes that have been caught in the LCR during 
any prior fall (i.e., potential residents) have a smaller 
chance of moving out of the LCR into the Colorado River 
than do adult chub that have never been caught in the LCR 
during the fall (Fig. 3B). Lastly, the monthly probabilities 
of moving between locations for subadults and juveniles 
are generally low (Fig. 3C). The lone exception is that a 
large proportion of juveniles emigrate from the LCR during 
the monsoon season. 

Survival and size transition rates 

Monthly survival rates are higher in the Colorado River 
than in the LCR for all size classes (Fig. 4A). In addition, 
survival generally increases for larger fish, with the excep- 
tion of survival for size class 5 fish in the LCR, which is 
estimated as lower than survival of size classes 3 and 4. 
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On the other hand, size transition rates are much higher 
in the LCR than in the Colorado River (Fig. 4B). To illus- 
trate the potential range of life histories that these sur- 
vival and size transition rates suggest, we consider two 

(A) Adult monthly movement rates into the LCR increase in 
the month prior to spawning, and a greater proportion of 
larger adults (size class 5) move in than smaller adults (size 
class 4). 



Colorado Monthly transition rate 



LCR 



Small adults during the month prior to spawning: 
0.31 (0.21 - 0.41) 



During rest of the year: 
0.004 (0.001 - 0.01) 



Large adults during the month prior to spawning: 
^^.^^^ 0-61 (0-49 -0-73) ^^^-v^,^^ 

^y^-X^^^ During rest of the year: '^^^^^'J^J-j^j^S^ 
0.04 (0.02-0.07) 

(B) Adults that have previously been caught in the LCR 
during the fall move out of the LCR at much lower rates 
than adults that have not been. 



Colorado Monthly transition rate 

Small adults with prior fall LCR capture 
All year: 0.01 (0.001 - 0.03) 



LCR 



Without prior fall LCR capture: 
2 months post spawning- 0. 17 (0.10 - 0.24) 
Rest of the year - 0* {estimated at zero) 

Large adults with prior fall LCR capture: 
0.06(0.01-0.11) ^ 

^ Without prior fall LCR capture: )yi^^^^ 
2 months post spawning- 0. 64 (0.2 — 0.93) 
Rest of the year - 0.40 (0.25 - 0.54) 



(C) Rates of movement in both directions are low for juveniles 
and sub-adults, with the exception that a large number of 
juveniles leave the LCR during summer monsoons. 



Colorado 



Monthly transition rate 

Juveniles 
0* (never observed during study) 



LCR 



During summer monsoons - 0. 61 (0.47 — 0.76) 
Rest of the year - 0.01 (0.003 - 0.03) 



Small sub-adults 
0.003 (0.002 - 0.005) 

0.01 (0.004- 0.03) 

Large sub-adults 
0.008 (0.005 - 0.011) 

0. 02 (0.007 - 0.03) 



Figure 3. Point estimates (and 95% confidence intervals) for monthly 
transition probabilities associated with three movement hypotheses. 



scenarios that bracket the potential variation: (1) life-long 
residency in the LCR starting in July of a chub's first year 
or (2) rearing in the Colorado River from July of the first 
year (and returning to the LCR only to spawn). Under 
the first scenario, chub quickly grow to adulthood and 
have short life expectancies as adults (Fig. 5). On other 
hand, chub that rear in the Colorado River take more 
than twice as long to reach adulthood and are -40% less 
likely to survive to be either small or large adults, but live 
~5 times longer, on average, as larger adults (assuming 
they reach adulthood). Under both scenarios, the 
expected number of spawns as a large adult (a proxy for 
fitness) is around 0.05-0.10. 

The best model estimates t as 0.18 (95% CI: 0.13- 
0.24), within the range we predicted a priori. Derived 
estimates of juvenile abundance in the Colorado River 
study site increased each year between July and September 
consistent with our juvenile movement hypothesis 
(Fig. 3C & Appendix S3). Derived estimates of other size 
classes and the overall adult population in the system 
were relatively stable over the study period, with the most 
precise estimate of total adult abundance coming in Sep- 
tember 2011 (11,000 95% CL 7000-16,000; Appendix S3) 



(A) 1 



0.95 



£ 0.9 



0.85 



• Little Colorado River 
^ Colorado River 



3 4 
Size class 



(B) 



• Little Colorado River 
A Colorado River 



0.20 
0.15 
0.10 
0.05 
i 0.00 

1,2 2,3 3,4 4,5 
Size class transition (from, to) 

Figure 4. Estimates of (A) monthly survival and (B) monthly size 
transition rates from best model (whiskers around dots indicate 95% 
confidence intervals adjusted by c). 
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Expected life history of a humpbaclc chub that spends its whole life in the LCR 



(A) Quick growth to adulthood 



Average number of 
months spent within 
size class 



Average inonthly 
survival rate 



0.9 



0.91 



0.96 



0.96 



(B) Short lives as large adults 

(average of 1-2 years as large adult) 



0.1 



■5 0.01 



^11 % chance of reaching 200 mm (from luly of year 0) 
>~5% chance of reaching 250 mm 



a. 0.001 




— I — 
20 



~~1 — 
30 



0 5 

Years as a large adult 



40 



Expected life history of a humpback chub that rears entirely in the Colorado River and spawns in the LCR 



(C) Slow growth to adulthood 



Average number of 
months spent within 
size class 




Average inonthly 
survival rate 



0.94 



0.97 



0.97 



0.98 

~4 % chance of reaching 200 mm (from luly of year 0) 
> ~2 % chance of reaching 250 mm 



(D) Longer lives as large adults 
(average of ~6 years as large adult 
if spawn 3 out of 5 years) 

D, 1 



If never spawned 



t 0.01 - 




Years as a large adult 



Figure 5. Comparison of two different strategies with similar overall fitness's. Chub that rear strictly in the LCR (A) mature quickly and (B) do not 
live long as adults. Chub that rear strictly in the Colorado River and then skip spawn (C) mature slowly, but (D) have much longer expected adult 
life spans. 



Discussion 

Our results are consistent with the existence of both 
shared breeding and skip-spawning partial migration in 
the chub population that spawns in the LCR (Table 2; 
Fig. 3). While skip spawning is known to occur in this 
population (Coggins et al. 2006b), we provide the first 
estimates of rates and show that larger adults in the Colo- 
rado River migrate into the LCR during the spawning 
season more frequently than smaller adults, consistent 
with theoretical expectations and observations in most 
other systems (j0rgensen et al. 2006). We also found that 
adults that were captured in the LCR during the 
nonspawning season are less likely to leave the LCR at 
any point in the rest of the study than adult chub that 
have not been captured there during the nonspawning 
season. This observation supports the hypothesis that 
some adult chub reside in the LCR year-round. 

Partial migration is believed to be relatively common, 
yet understudied, and we are only beginning to under- 
stand how humans are affecting these systems (Chapman 
et al. 2011). Here, we have shown that there can be large 
differences in vital rates within different parts of the river 
system that includes both unregulated and regulated habi- 
tats, and these differences in vital rates can lead to very 
different life histories for individuals that adopt a strictly 



migratory or strictly breeding ground resident strategy 
(Fig. 5). Despite the large differences in growth and 
survival between the Colorado River and the LCR, the fit- 
ness of the resident and migrant strategies may be 
roughly equal. For example, if we chose the number of 
spawns as a large adult as a proxy for the fitness of a 
given strategy, then expected fitness for a juvenile in July 
does not differ significantly based on whether it resides in 
the LCR for the rest of its life (0.05-0.1 expected spawns 
based on the probability of reaching large adulthood and 
life expectancy as a large adult: Fig 5A-B) or moves into 
the Colorado River and rears until adopting a migratory 
strategy as an adult (-0.07 expected spawns based on 0.02 
probability of reaching large adulthood, the expected life 
span of 6 years for a skip-spawning adult, and spawning 
rate of 0.6: Fig. 5C-D). On the other hand, the genera- 
tion time of residents is roughly half that of migrants, 
which may mean both quicker declines and recoveries of 
residents relative to migrants. 

While intriguing, these calculations are nonetheless 
based on a short period of time (3.5 years for a species 
that can live to be 30 -I- years), ignore substantial uncer- 
tainty, and occurred during a period when conditions in 
the Colorado River were favorable. For example, average 
water temperatures in the regulated Colorado River were 
~1°C colder during the 1980s through early 2000s as 
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compared to our study period (and presumably chub 
growth rates were hkewise lower). Furthermore, salmonid 
abundances near the confluence were relatively low dur- 
ing the study period compared with estimates from other 
times in the last few decades (suggesting that HBC sur- 
vival rates in the Colorado River may have been substan- 
tially lower during periods of time in the past; Yard et al. 
2011; Korman et al. 2012). It is impossible to determine 
how present variation in chub growth and survival in the 
Colorado River compares to the conditions over evolu- 
tionary time periods given changes in many factors 
expected to affect these rates including introduction of 
nonnative fish species, extirpation of native fish species, 
shifts in the food base and massive changes in the physi- 
cal template (especially temperature, flow and turbidity). 

While the regulated Colorado River has changed drasti- 
cally over the last century, the LCR is also modified. 
From the perspective of chub, the greatest modification 
has been the addition of nonnative species that prey on 
humpback chub (Marsh and Douglas 1997). Determining 
abundances of many nonnative species in the LCR is 
difficult because they are not easily captured using hoop- 
nets. The impacts of these nonnative species on chub sur- 
vival are unknown, and it has been suggested that they 
could be significant, particularly with respect to smaller 
size classes. While we had no a priori hypotheses concern- 
ing differences in survival rates between the LCR and the 
Colorado River, we were surprised that monthly survival 
rates were consistently lower among all size classes in the 
LCR. In particular, the apparent decline in estimates of 
survival between size class 4 and size class 5 fish in the 
LCR was unexpected and will be a focus of future 
research. One a posteriori hypothesis is that larger adults 
in the LCR are experiencing increased mortality associ- 
ated with the physiological costs of spawning, while a sec- 
ond a posteriori hypothesis it that there is long-term trap 
response occurring in chub that have been caught repeat- 
edly. Future analyses using data from arrays of PIT tag 
antennas spanning the width of the LCR may help test 
the second hypothesis. 

While we find support for adult movement consistent 
with the residents and skip-spawning migrants present in 
theoretical models of partial migration, we also note low 
rates of movements between the LCR and Colorado River 
that suggest additional complexity. Thus, while some por- 
tions of the chub population are likely to show life histo- 
ries similar to the extremes of life-long residency and 
rearing in the Colorado River and only returning to the 
LCR for spawning Olustrated in Fig. 5, other portions of 
the population are likely to show intermediate growth rates 
(and survival). Likewise, while rates of movement between 
the LCR and the Colorado River for subadults and 
juveniles are generally low, they stiU suggest some move- 



ment. These low monthly rates compound over the course 
of a year and can lead to substantial annual fluxes of indi- 
viduals. For example, when population sizes are in the 
thousands, a monthly movement rate around 0.01 could 
lead to a flux of hundreds of individuals over the course of 
a year. 

The one exception to low monthly rates of movement 
prior to adulthood is the high probability of juveniles 
moving from the LCR to the Colorado River during the 
monsoon season. The effects of these movements on 
abundances of juveniles in the Colorado River study site 
are particularly clear in 2011 and 2012 when estimates of 
juvenile abundances increase by -4000 between July and 
September of each year (Appendix S3). The increase in 
abundances between July and September of 2010 (and 
between August and September 2009), in contrast, is 
much smaller, suggesting either that the export rate is not 
constant or that there was a smaller population of juve- 
niles in the LCR to begin with in July. Interestingly, there 
was substantially more flooding during the monsoon sea- 
son in 2011 and 2012 than in 2009 and 2010, suggesting 
a potential link between export and LCR flood frequency. 
Unfortunately, we are not able to estimate juvenile abun- 
dances in the LCR during July, and the relatively small 
number of marks put into the LCR during July does not 
allow us to estimate year-specific rates of outmigration 
over the study period. Ongoing research is increasing 
effort in July to obtain abundance estimates and allow for 
yearly estimates of outmigration rates. 

An unresolved question in our river system regards the 
degree to which the resident versus the migrant popula- 
tions drive overall dynamics of the partial migratory sys- 
tem. Based on our estimated rates, we expect that the two 
strategies would lead to similar reproductive results; how- 
ever, we might also expect that the resident strategy would 
outperform the migrant strategy when the Colorado River 
is less hospitable (e.g., higher salmonid abundances or 
cooler water temperatures). As LCR residents would have 
a quicker generation time, they should also increase (or 
decrease) more quickly. On the other hand, 82% of the 
adult population (based on abundances reported in 
Appendix S3 and our estimate of t) currently resides in 
the Colorado River during the nonspawning season, and 
vital rates under the environmental conditions experi- 
enced in 2009-2012 suggest that juveniles can successfully 
rear to adulthood in the Colorado River mainstem 
(Fig. 5). Survival was once thought to be uncommonly 
rare in the Colorado River because of the seasonally 
constant, low temperatures (Clarkson and Childs 2000). 

One explanation is that under current environmental 
conditions, both strategies have roughly similar fitness's 
and are contributing to a fairly stable or slightly increas- 
ing adult population; however, when environmental 
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conditions in the Colorado River decline, the migrating 
population in the Colorado River may become a popula- 
tion sink (sensu PuUiam 1988). Given the limited carrying 
capacity in the LCR (Meretsky et al. 2000; Benenati et al. 
2002), this suggests that population targets in the system 
can only be met by managing the Colorado River to 
ensure relatively beneficial environmental conditions (e.g., 
lower salmonid abundances or warmer water tempera- 
tures). Moreover, as system dynamics are relatively slow, 
conditions in the Colorado River will likely have to per- 
sist for long periods of time to change adult population 
sizes. Vital rates and abundances of smaller size classes 
will provide a leading indicator of future adult population 
size. However, given the long periods of time spent in 
these vulnerable life stages, a year or two of beneficial 
conditions can be easily negated by a change in these con- 
ditions. Future work will focus on exploring annual varia- 
tion in vital rates, particularly as salmonid abundances in 
the Colorado River may be increasing. 
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